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We investigate the QCD phase diagram by using the strong-coupling expansion of the lattice QCD with one species of staggered 
fermion and the Polyakov loop effective action at finite temperature (T) and quark chemical potential (//). We derive an analytic 
expression of effective potential f e s including both the chiral (17^(1)) and the deconfinement (Z^ ) dynamics with finite coupling 
effects in the mean-field approximation. The Polyakov loop increasing rate {d( p /dT) is found to have two peaks as a function of 
T for small quark masses. One of them is the chiral-induced peak associated with the rapid decrease of the chiral condensate. The 
temperature of the other peak is almost independent of the quark mass or chemical potential, and this peak is interpreted as the 
Z Nc -induced peak. 



Keywords: Lattice QCD, Extreme QCD, Strong-coupling, Deconfinement, Chiral Symmetry 



D ■ The phase diagram for chiral and deconfinement transitions 
, jn Quantum Chromodynamics (QCD) at finite temperature (T) 
and/or quark chemical potential (//) is one of the most fascinat- 
t-H ing subjects in the current high energy and nuclear physics. To 
k* - investigate properties of QCD at high T is one of the physics 
^| goals in the LHC-ALICE experiments. At /i — 0, first princi- 
pie investigations based on lattice Monte-Carlo (LQCD-MC) 
^—1 simulations predict pseudo-critical temperature, T c 145 - 
•195 MeV, for the chiral phase transition U|. In the finite /j 
region, the LQCD-MCs suffer from the notorious sign problem 
of the quark determinant, and have not provided reliable results 
in the large chemical potential region, fi/T > 1. 
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The strong-coupling (l/g 2 ) expansion in the lattice QCD 
(SC-LQCD) has been successful since the beginning of the 
lattice gauge theory, and would provide an alternative lat- 
tice framework to study the QCD phase diagram including 
finite // region. In pure Yang-Mills theory, the string ten- 
sion in the strong-coupling limit gives the area law [2], and 
the LQCD-MC yfl smoothly connects the strong-coupling re- 
sult 14:!] to the scaling region. In the pure Yang-Mills theory 
at finite T, we can describe the Zn c deconfinement transition 
based on the effective potential for the Polyakov-loop {£ p ) in 
the leading order of the strong-coupling expansion with the 
Haar measure effects Jl, and higher order corrections have 
been investigated recently [H. For the SC-LQCD including 
fermions, many theoretical knowledge have been accumulated 
so far 0, S H OH O, [H [H Q El, and the chiral phase 
transition in the T - u plane have been well investigated in the 
strong-coupling limit |3 17, 18, 19]. It is remarkable that the 
coupling of the chiral condensate cr and the Polyakov-loop £ p 



was extracted in the strong-coupling limit [20[ 21], and led to 
the invention of the Nambu-Jona-Lasino model with Polyakov- 
loops (PNJL model) 12211 . Recently, the finite lattice couping 
(J3 = 2N c /g 2 ) effects are incorporated, and is found to give rise 
to modifications of quark mass and chemical potential. The 
evolution of the chiral phase transition with increasing f3 has 
been interpreted via these modifications J23, 24]. As will be 



shown later, this development opens a possibility to investigate 
the finite coupling evolution of Zjy e deconfinement dynamics in 
addition to the chiral dynamics in the whole T — n plane 12511 . 

One of the interesting observations in the LQCD-MC is 
that peak positions of chiral and Polyakov loop susceptibili- 
ties (Xo-,c) are close to each other, and the small separation of 
them could be explained as a consequence of the broad ana- 
lytic behavior of the crossovers [Q]]. It would be meaningful to 
ask ourselves whether the peak of Xi is induced by the chiral 
crossover. Otherwise, does the Zjq c dynamics (accidentally or 
inevitably) leads to the peak oixi , near the chiral crossover? To 
shed light on this problem, it is a good strategy to investigate the 
finite /u cases. In models such as PNJL model combined with 
the statistical model 12611 or Polyakov-Quark Meson model with 
the functional renormalization group evolution 02711 . two transi- 
tions almost coincide. In the PNJL with a certain fit parameter 
set, the first-order chiral phase transition with a small jump of 
a small value of t„ can be realized in the low T and large yt 
region SHIES (Sh- 



in this Letter, we investigate the chiral and Z^ c deconfine- 
ment dynamics by using the SC-LQCD with the Polyakov loop 
effects, abbreviated as P-SC-LQCD. The P-SC-LQCD is di- 
rectly based on the lattice QCD, and it does not contain any 
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additional parameters than those in QCD. The lattice coupling 
P = 2N c /g 2 in the plaquette action is a unique parameter of 
the lattice QCD in the chiral limit. To investigate the chi- 
ral and deconfinement dynamics simultaneously, we con- 
sider the effective action with leading [0(1 /g )] and next-to- 
leading order [NLO, 0{1/ g 2 )] effects of the strong coupling 
expansion in the fermionic sector, and the leading order con- 
tributions to the Polyakov-loop [0(1 / g 2N ')] in the pure gluonic 
sector. The present framework aims at developing previous SC- 
LQCD studies for the chiral dynamics |Q~0, El GS H 124 ] 
to include the Z^ c deconfinement dynamics, and leads to an ex- 
tended version of the IK-GO model (2(1 21 1 to include finite 
beta effects for the quark sector. This framework allows us to 
investigate the beta evolution of the interplay between chiral 
and Zjv c dynamics consistently. This is the advantage of using 
the P-SC-LQCD over effective models. 

We briefly overview the derivation of the effective potential 
in P-SC-LQCD. Details are shown in our previous papers Q23 



X + J + k 



3211 . We start from the lattice QCD partition function with one 
species of staggered fermion (x) with a current quark mass (mo) 
in the lattice unit a — 1 , 



S F - S G - m Q ^Xx 



Zlqcd = J £>\X>X, U v ] exp 

Sf = ^ ^[r]v,xXxU v , x )c x+v - T]-*(h.c.)] 

v,x 
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where U v , x e SU(N C ) and U P= ^ X = ^[U^Uy^ul^ulj 
represent the link- variable and plaquette, and the staggered sign 
factor j] YX = exp(jU(5 v o)(-l) ;to+ '" +x '" 1 contains the the lattice 
chemical potential fi. The transformation Xx — > e' ee *Xx with 
£ x - (—l) Xo+ "' +Xd turns out to be a L^l)^ chiral transforma- 
tion UGll, which leaves the staggered action S f invariant in 
the chiral limit (mo — > 0). The U x (l) chiral symmetry would be 
enhanced to SU(N f = 4) in the continuum limit JTfl, 0, H ■ 
The plaquette action Sc is invariant under the global trans- 
formation U ViX — > ClU ViX , where Q. is the element of the cen- 
ter of the SU(N C ) gauge group, Z^ c . The chiral condensate 
(cr ~ (xx)) and the Polyakov-loop (L p = tr e (n T U 0iST )/N c ) are 
the order parameters of the chiral and Z^ c symmetries, respec- 
tively, and they take finite values via the spontaneous breaking 
of these symmetries. Both chiral and Z^ c symmetries are ex- 
plicitly broken for a finite quark mass mo, and the transitions 
could be replaced with the crossovers. We concentrate on the 
color S\J(N C = 3) in the 3 + 1 dimension (d = 3) in the later 
discussion. 

The effective action of hadronic composites is obtained by 
the Taylor expansion in fi and integrating out spatial link vari- 
ables in the finite T treatment of P-SC-LQCD. We include the 
leading and NLO terms in the fermionic sector, ll23ll . and 
we adopt the leading order Polyakov-loop effective action Sg 
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Figure 1: Effective action terms in the strong-coupling limit and l/g 2 correc- 
tions. Open circles, Filled circles, and arrows show x,X, and U r , respectively. 



in the pure gluonic sector 13511 . 
Zlqcd - J !D[x,X, Uol exp (-S^° - S™) , (4) 
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S™[L p ,l p ] = -AT^ J_p 1/r J][L p , x L p ^j + h.c] , (6) 

° c M 

where the coupling /3 T . S and composites are defined as (j3 T , /3 S ) = 
(/3d/2N 3 c , /3d(d~ l)/16iV|), M x = XxXx, and (V», V;Qi)) = 
(e^XxU^xXx+fr e^Xx+oU^Xx)- We consider only the leading 
order terms in the l/d expansion |@], which corresponds to the 
minimum quark number diagrams for a given plaquette config- 
uration as shown in Fig.Q] «S^° contains the four- Fermi inter- 
action, the temporal components of the vector interaction, and 
eight-Fermi interaction, and shows the nearest-neighbor 
interaction of the Polyakov loops. We can reduce the effective 
action S^f° into the bi-linear form in staggered fermions by in- 
troducing several auxiliary fields (cr, u> T , and (p TtS ) summarized 
in Table Q] 
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G^(m q ,P) =ih q 5 xy + -f(e*Uo, x 6 x+h - e ,«)., 6y ) , (8) 

where N T (N s ) represents the temporal (spatial) lattice size, 
and all auxiliary fields are assumed to be constant and static. 
The spontaneous breaking of the chiral symmetry with NLO 
effects results in the dynamical shifts of quark mass iii q = 
[mo + (d/2N c + 2fis<p s )o~] /Z x , chemical potential jj. = fi - 
\ogZ x , and the quark wave function renormalization factor, Z x , 
where Z x = VZ+Z and Z ± = 1 + fi T (<p T ± <jd t ). Such a prop- 
erty allows us to utilize the technique developed in the strong 
coupling limit, and from a technical point of view, it is essential 
to derive the effective potential including both the chiral and 
deconfinement dynamics. 



Table 1: The auxiliary fields and their stationary values. Here, tpo = N c -Z x m c/ + 
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Figure 2: Schematic figure of the Polyakov-loop construction in the strong- 
coupling expansion. 



Now we perform the Gaussian integral over the staggered 
quark fields (x,x)- Here we work in the static and diagonalized 
gauge (called the Polyakov gauge) for temporal link variables 
with respect for the periodicity 111 ill . 



'Wo,* = f[ U °* = diag{ e ' fl ' • • • , 



(9) 



Owing to the static property of the auxiliary fields and the tem- 
poral link variable in the Polyakov gauge, the quark determi- 
nant is factorized in terms of the frequency modes, and evalu- 
ated by the Matsubara method (see for example the appendix 
in Ref. Il8ll ). The remarkable point is that the resultant expres- 
sion can be expressed in terms of the Polyakov-loop variables 
(L P ,L P ), 

I T>\x,x\ exp -Z x ^ XxG:,l {m q ,p)Xy 

J XV 

~J[e N ^' T K q (T, M )'R q (T, t x)] 



7 N c N T N c ! 



K q (T,n)= 1+e 



-N c (.E q -p)IT 



R- q (T,n)= l+ e - N ^ +il) ' T 
+ N C (L 



-{E q +fi)IT r -2{E q +ji)IT\ 
p,x c ' U PJ& ) • 



(10) 

(11) 
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where E(m q (cr)) - sinlT 1 [m 9 (<x)] corresponds to the quark ex- 
citation energy. In that expression, L p couples to a Boltzmann 
factor e~( £ « _ A)/ :r j an( j determines how quarks thermally excites. 
In the confined phase (L p ~ 0), color-singlet states dominate, 
while quarks can excite in the deconfined phase (L p + 0). Equa- 
tions ( fTTT ) and ( fT2l give a natural coupling manner between L p 
and cr. This point has been pointed out in the strong-coupling 



limit 1201 12 111, and utilized in the PNJL model |22j|. In the cur- 
rent formulation, the Boltzmann factor includes the NLO ef- 
fects in (m q ,p), and the coupling manner between L p and cr is 
modified by the finite coupling effects. 

Finally, we evaluate the temporal link integral and obtain the 
effective potential. In the Polyakov gauge, the Haar measure 
becomes a Van der Monde determinant over the color space, 
and can be rewritten by using the (reduced) Polyakov-loop 
(L P ,, -» Z^i e w ^IN c ), 

j fl"l/ ,x = 27 J d[L p x , L PtX ]Mliaar(L p , L p ) , (13) 



M H ^(L P ,L P ) = 1 - 6L P , X L AX - 3(L p , x L p , x y + 4{L% + L%) 



(14) 



While it is possible to perform this integral exactly, we here 
adopt a simpler prescription; we replace Polyakov-loops in the 
integrand with its constant mean-field value, (L PiX ,L PjX ) — > 
(t p ,{p), and search for the stationary values of (t p ,( p ). This 
treatment gives the effective potential in a similar expression to 
that used in the PNJL model, and useful for the comparison. 
Thus, we obtain the effective potential as a function of the aux- 
iliary fields <Tj = (cr, (p T S , u) T , l p , lp), temperature T, and quark 
chemical potential /i in the mean-field approximation, 

^ eff (<D; T, n) = -(T log ZlqcdWs = ?? ff + ^ , (15) 

„2 



- I 4r + ^ s y + P -f- + f - col) - N c logZ, , 



eff \4N C 

- N c E q - TQogKiiT.n) + log^(r,Ai)) 



(16) 



T^ x - -2TdN 2 {^y~{ p { p - T log Mu^dp, € p ) , (17) 



The equilibrium is determined by imposing stationary condi- 
tions on the effective potential, dT~ e ff / <3<t> = 0, which lead to the 
relations summarized in the third column of Table [TJ T* s is re- 
sponsible for the chiral-dynamics, and T^ 1 originates from the 
plaquette action and governs the Zn c dynamics. These ingredi- 
ents communicate with each other through the quark determi- 
nant effects i.e. % q {T,n) and ^(T,/i) in Eq. (O. 

The I pip term in Eq. ( TTTb gives large finite T effects to Zjq c 
dynamics at finite f3, and vanishes in the strong-coupling limit. 
In the previous work in the strong coupling limit, this quadratic 
term is fixed to a constant to be consistent with the empirical 
value of the string tension 12 ill . The finite coupling property of 
the current formulation allows us to investigate the f3 evolution 
of the Iptp term, i.e. finite T effects of Zjy c dynamics, without 
introducing additional parameters. 

We shall show the first P-SC-LQCD results including the 
"chiral and deconfinement" dynamics, "finite /i effects", and 
"finite coupling effects" simultaneously. As shown in our re- 
cent work [|32j], the critical temperature at zero chemical poten- 
tial becomes closer to the LQCD-MC results in the coupling 
region /? < 4. Therefore, we focus our attention to the results 
at f3 — 4. In the last part, we discuss the stability of our main 
conclusion for variations of f3. 
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Figure 3: The phase boundary for the chiral transition with the second peak of 
dCp/dT at /? = 4 in the chiral limit. The "P" and "Q" correspond to those in 
Fig.|4]and in the upper panel of Fig. [6] 



In Fig. [3j the phase diagram at (J3, mo) = (4, 0) is shown 
in the lattice unit. We find the first- (solid blue) and second- 
order (dashed blue) chiral transition lines separated by the (tri- 
)critical point (CP) at (// C p,7cp) = (0.58,0.19). The result is 
qualitatively consistent with the previous SC-LQCD with NLO 
effects II23I1 . In Fig. [4] we show the T dependence of the chi- 
ral condensate cr and the Polyakov loop t p . The upper (lower) 
panel displays the results for jj — (0.5), i.e. on the T axis 
(dash-dotted line) of the phase diagram in Fig. [3] We note 
jj = 0.5 < /icp, and the following results would not be con- 
taminated by the fluctuation around the CP. 

We find two peaks in dt p /dT. One peak appears at the chiral 
phase transition. For both fi = and 0.5 cases, the Polyakov 
loop C p shows a rapid increase at the chiral transition tempera- 
ture (cr — > 0). The strong correlation of the Polyakov loop and 
cr is found at any point on the chiral transition boundary. This 
correlation can be seen more clearly in the derivative dt p /dT as 
shown in Fig. [5] We find a sharp peak in the vicinity of the chi- 
ral phase transition, T c ^ i= q(J3 = 4) 0.486. The almost simulta- 
neous observation of the chiral transition and the rapid change 
of t p would be consistent with the LQCD-MC results d. Via 
the cr-lp coupling in the Boltzmann factor terms in Eq. (fTTT i and 
dl2l . the chiral transition would induce the rapid change of t p . 
Thus we regard this peak as the chiral-induced peak. 

We also find the second peak above the chiral transition tem- 
perature in dtp/dT. At fi = 0, a small enhancement is seen 
at around T - 0.52. At finite ji, the enhancement of dt p /dT 
at the second peak becomes significant. A similar double-peak 
structure has been reported in the model studies based on PNJL 
model l29ll . In Fig. |6l we show dt p /dT along the line with 
p = 0.5 (dash-dotted line in Fig. [3}. In the chiral limit (mo — » 0, 
upper panel), we find two peaks "P" and "Q". Here, "P" and 
"Q" in Fig. [6] correspond to those in Figs. [3]and|4] The peak 
"P" locates at the chiral phase transition, and is clearly inter- 
preted as the chiral-induced peak. As indicated by the red line 
in Fig. |3] the similar peak to "Q" is observed in the whole range 
of yu. The strength of this peak becomes weaker for smaller fi, 
which is expressed by the dotted line in Fig. [3] 

The peak "Q" can be understood as a signal of the Z^- 
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Figure 4: The T dependence of the Polyakov loop and the chiral condensate 
for n = (upper) and 0.5 (lower) at (fi, mo) = (4, 0). In the lower panel, the "P" 
and "Q" correspond to those in Fig.|3]and the upper panel of Fig. [6] 



induced crossover: Zjy c is the symmetry in the pure gluonic sec- 
tor, and it becomes exact in the heavy quark mass limit. We also 
expect weak fi dependence of the Z^ c deconfinement transition, 
since does not directly depend on quark chemical potential 
yU. In the middle and lower panels of Fig. [6] we show dtp/dT for 
mo = 0.03 and mo = 1, respectively. For mo — 0.03, we find two 
peaks. The temperature of the first peak "P" is slightly shifted 
upward and becomes closer to "Q" (Tq) with increasing mo, 
while Tq stays almost constant. For larger masses, mo > 0.05, 
the two peaks merges to a single peak, as shown in the lower 
panel of Fig.|6]formo = 1. This single peak grows with increas- 
ing mo, and its temperature is nearly mo independent, T ~ 0.52, 
which is close to Tq at smaller quark masses. The -induced 
nature of "Q" and the merged single peak are confirmed by the 
weak dependence of Tq on fi and mo as found in Figs. |3]and|6] 
respectively. It is interesting to find that the Z\ c nature survives 
in the chiral limit or small mass region, and can be observed as 
a peak in dt p /dT. 

There are several comments in order, (a) In cold dense matter 
in Fig. [3] the chiral symmetry is restored and Polyakov loop is 
suppressed. These features may be similar to those of quarky- 
onic matter ll30h . (b) We find that Tq is larger than the chiral 
phase transition temperature for small quark masses. In the chi- 
ral limit, the Polyakov loop f p decouples from the chiral dy- 
namics at T — Tq since <T-t p couplings in Eq. <TT3T > vanish due 
to the exact chiral restoration cr — 0. (c) We note that the chiral 
condensate decreasing rate, —dcr/dT, has a small peak at the 
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Figure 5: The T dependence of dC p /dT in the chiral limit with zero chemical 
potential for fl = 4. The values are shown in the lattice unit a = 1. The dashed 
blue line represents the critical temperature T c fl= o for the second-order chiral 
transition. 



Zty, -deconfinement crossover for mo — 1 as shown in the lower 
panel of Fig. [6] This peak is interpreted as a -induced peak. 
In this case, however, another peak which would stem from the 
chiral symmetry is completely overwhelmed due to the large 
quark mass, and the double-peak structure does not appear. 

Finally, we discuss the /3 dependence of the two peak struc- 
ture of dtpldT at finite fi. In Fig. |7j we show d{ p /dT for 
several (3 at (mo,fi) = (0.01,0.5). As indicated from this fig- 
ure, the two peaks are found at least in the range 2 < (3 < 6. 
For each /?, we define the pseudo-critical temperatures for the 
chiral-induced (T^ _ 5 (J3)) and Z^ c -induced (7^ _ 5 (y6)) decon- 
finement crossovers as the first and second peaks of d£ p /dT, 
respectively. Both of them are decreasing functions of /3. For 



r^ =O5 (j0), the current results are close to the pseudo-critical 
temperatures for the chiral phase transition in our previous 
works J23, 24]. This again indicates the chiral induced nature 



of the first peak "P". 

We find that the dependence of T^_ 05 (fi) is larger than 

that of T^^sfJS), and the separation between two peaks tends 
to be narrower with increasing /? when chemical potential is 
fixed, fi = 0.5. Whereas the two peaks tend to be more sepa- 
rated for larger chemical potential for a fixed value of /3. For 
example, two peak separation at fi = 0.55 becomes 1.4 times 
larger than that at fi — 0.5 for j3 — 6. This is because the chiral- 
induced transition temperature decreases in the large fi region 
of the T - fi plane, while there is no direct fi dependence in the 
Zff c dynamics. 

In summary, we have investigated the chiral and deconfine- 
ment crossovers at finite temperature T and quark chemical po- 
tential fi based on the strong-coupling expansion in the lattice 
QCD with one species of staggered fermion. We have consid- 
ered the leading and NLO effects in the strong-coupling ex- 
pansion for fermionic sector, and the leading order Polyakov- 
loop effective action terms for the pure Yang-Mills sector. The 
Z\ c deconfinement dynamics has been incorporated to the Haar 
measure, where the Polyakov-loop is replaced with its constant 
mean-field value. 

We have found double-peak structure in the Polyakov loop 
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Figure 6: T dependence of d( p /dT 3.1)1 = 0.5 for /? = 4 in the lattice unit. The 
current quark masses (mo) are 0.0 (upper), 0.03 (middle), and 1 (bottom). The 
dashed lines in the middle and lower panels represent -(iV c )~' dtr I dT . 



increasing rate dt p /dT as a function of T for small quark 
masses and large fi. The first peak is induced by the chiral 
transition. This is because the Polyakov loop t p becomes sen- 
sitive to the chiral dynamics through the coupling to the chiral 
condensate <x. For the larger quark mass mo, the first peak is 
overwhelmed by the second peak, whose position is almost in- 
dependent of mo. This indicates that the second peak would 
attribute to the remnant Z^ c dynamics, which is less affected by 
fi than the chiral phase transition line due to the lack of a direct 
fi dependence in the Haar measure treatment for Z^ c . Hence the 
double peaks, i.e. the chiral and Z^ c induced peaks, come out 
in large fi region. 

As future perspectives, we should evaluate the chiral and 
Polyakov-loop susceptibilities, next-to-next-to-leading order 
(NNLO) terms of strong coupling expansion for the chiral dy- 
namics 12411 . higher order corrections of the Polyakov-loop ef- 
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Figure 7: The T dependence of d{ p /dT for several f3 at (mo.fi) = (0.01, 0.5). 



fective action JfJ, and higher order terms of the 1 Id expansion. 
The exact evaluation in each order of the strong-coupling ex- 
pansion is also expected by the Monomer-Dimer-Polymer for- 
mulation 1 15 , . These corrections include couplings between 
the Polyakov loop and the chiral sector in the effective action 
level. Taking account of the resultant entanglement effects of 
the chiral and Z^ c dynamics, the appearance of chiral and Ztf c 
induced peaks must be investigated in future. 
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